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Abstract 

An exact method is suggested to treat the nonlinear self-interactions (NLSI) in the 
relativistic Hartree-Fock (RHF) approach for nuclear systems. We consider here the 
NLSI constructed from the relativistic scalar nucleon densities and including products 
of six and eight fermion fields. This type of NLSI corresponds to the zero range limit 
of the standard cubic and quartic self-interactions of the scalar field. The method to 
treat the NLSI uses the Fierz transformation, which enables one to express the exchange 
(Fock) components in terms of the direct (Hartree) ones. The method is applied to 
nuclear matter and finite nuclei. It is shown that, in the RHF formalism, the NLSI, 
which are explicitly isovector- independent, generate scalar, vector and tensor nucleon 
self-energies strongly density-dependent. This strong isovector structure of the self- 
energies is due to the exchange terms of the RHF method. Calculations are carried out 
with a parametrization containing five free parameters. The model allows a description 
of both types of systems compatible with experimental data. 

PACS number(s): 21.30.Fe, 21.10.-k, 21.60.Jz, 21.65,+f. 

Keywords: relativistic Hartree-Fock approach, Fierz transformation. 
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Short title: Relativistic Hartree-Fock description of nuclear systems. 



1 Introduction 



The Relativistic Hartree-Fock (RHF) approach for finite nuclei has been developed in Refs. 
[1-8] (see also references therein) for the so-called linear models, which are characterized by 
linear field equations. In that approach, the contribution of the pion degrees of freedom 
(single pion exchange) is taken into account explicitly. In Ref. [2] it is shown, in particular, 
that the pseudovector coupling of pions to nucleons is more preferable than the pseudoscalar 
one in the nuclear structure context. In Ref. [9], certain important features concerning the 
spin-orbit interation are directly related to the pion contribution. Thus, the incorporation of 
pions into the model represents one of the main advantages of the RHF method in comparison 
to the relativistic Hartree approach. However, in the papers mentioned above the nonlinear 
self-interactions (NLSI) of the mesonic fields have not been taken into account. Including 
nonlinear self-interaction terms corresponds to one of the possibilities to account for the 
three- and four-body forces in the nuclear structure calculations. Different types of self- 
interaction Lagrangians have been considered in literature up to now. Initially, the scalar 
field self-interactions have been introduced involving a 3 - and <j 4 -terms, where a corresponds 
to the nuclear scalar field [11]. This type of self-interactions has been shown to play a very 
essential role in the relativistic Hartree calculations [12,13] to get, for example, the proper 
value of the compressibility modulus {K). One of the main problems of the RHF approach 
is that it brings about a high value of K. Thus, the inclusion of the NLSI terms in the 
RHF famework could solve this problem too. To work out the exact RHF equations with 
the NLSI is a complicated task not solved up to now. In Ref. [14], an approximate method 
to take into account self-interactions of the a 3 - and cr 4 -type in the RHF procedure has been 
suggested (see Ref. [15] for other type of NLSI). This method involves a simple idea based 
on the inclusion of the nonlinear terms appearing in the equation for the a field, together 
with the scalar meson mass m CT , into a scalar meson effective mass m*, which replaces m a 
in the corresponding meson propagator in the nuclear medium. Let us mention also that in 
Ref. [10] the authors present a procedure which allows one to evaluate the contribution of 
the Fock terms in a truncation scheme, including self-interactions of the scalar meson field, 
and illustrate their method in the case of infinite nuclear matter. 

In the present paper, we study the properties of a Lagrangian including, beside the 
exchange of cr, uj, it and p mesons between nucleons, the NLSI of the scalar field in the zero 
range limit (ZRL), in the framework of the RHF approximation. The ZRL allows one to 
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express the exchange (Fock) terms for the NLSI explicitly via the direct (Hartree) terms in 
an exact way. 

The paper is arranged as follows. In Sect. II, the general formalism is presented 
extensively: the Hartree-Fock contributions of the NLSI are calculated (in the ZRL without 
any approximation), both to the total energy of the system and to the nucleon self-energies, 
the expressions obtained being valid for nuclear matter (NM) and finite nuclei. For both 
cases, the results are shown in Sect. III. Finally, the conclusions are drawn in Sect. IV. 



The effective Lagrangian of the model considered in the present paper involves interactions 
of the nucleons via the exchange of mesons with the space-time transformation properties of 
the scalar a, vector (both isoscalar u, and isovector p), and pseudoscalar 7r fields. It contains 
a "linear" part, which generates linear field equations identical to the corresponding linear 
part of Ref. [14], and it will not be reproduced here. Let us mention only that the pion field 
is coupled to nucleons through a pseudovector coupling [2]. The current Lagrangian involves, 
however, a self-interaction part that is somewhat different from its homologous part of Ref. 
[14]. It takes the form: 



where ip is the fermion field operator, b and c are the coupling constants of the nonlinear 
terms, g„ and m a are, respectively, the scalar meson coupling constant and mass. We shall 
utilize also the dimensionless coupling constants b = and c = M being the nucleon 
mass. 

Notice that the NLSI given in Eq. (1) coincide with the conventional scalar field NLSI 
used in Ref. [14] in the ZRL of the scalar field, i.e., when the term with ml in the equation 
of motion of the scalar field dominates over the Laplacian and nonlinear terms. This is what 
happens for nuclear systems at small densities and with a smooth surface. Although, strictly 
speaking, real finite nuclei do not satisfy these conditions, for values of m a ~ 500 MeV, we 
can consider that this approximation in effective Lagrangians is acceptable. The parameters 
b, c, b and c have the same units as in Ref. [14] . 

Let us mention also that the NLSI given by Eq. (1) have just the same structure as 



2 General formalism 



A. The effective Lagrangian 
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the respective components of the Lagrangian appearing in the point coupling model and 
containing higher order terms in the fermion fields [17,18]. The inclusion of this kind of 
terms in the Lagrangian can be also justified by the necesity to introduce an extra density 
dependence in the Dirac-Brueckner-Hartree-Fock calculations to allow a simultaneous fit to 
the NN phase shifts and the nuclear matter equilibrium point [19]. 

B. The Hartree-Fock equations 

To get the RHF equations corresponding to our effective Lagrangian, we closely fol- 
low the Refs. [3,14], restricting ourselves to the static approximation for the meson fields. 
The nucleon field ip is expanded into a complete set of stationary single-particle spinors 
{ip a (x)e~ tEat } and we consider the tree approximation. 

The Dirac equation including exchange terms can be obtained by minimizing the total 
energy of the system, which is given by the expectation value of the total Hamiltonian in 
the space of Slater determinants il^a+lO >, where a+ is a creation operator for a nucleon in 
the state a. 

B. 1 Contributions of the NLSI terms to the total energy of the system 

The contribution to the energy of the linear part of the Lagrangian, of finite range, 
is taken directly from Ref. [3]. It includes, of course, the exchange terms. Here, we shall 
concentrate in the technique of calculating direct (Hartree) and exchange (Fock) components 
of the contribution to the energy of the self-interaction term Usi given by Eq. (1). 

In doing this, we obtain for the mean value of the cubic component of Usi 



< 0|(^) 3 |0 > = (0|£ £ ^(l)^(2)^ fc (3)^K3)^(2)^(l)a t +ata+a fc /a J va i ,|0) 

ijk i'j'k' 

= £^(1)^(1)^(2)^(2)^(3)^(3) 

ijk 

-3 £^(1)^(1)^(2)^(2)^(3)^(3) (2) 

ijk 

+2 £^(1)^(1)^(2)^(2)^(3)^(3), 

ijk 

and for the mean value of the quartic component 

<0|(^) 4 |0> =(0|£ £ ^(l)^(2)^ fc (3)^ m (4)^ m /(4)^(3)^-/(2)#(l) 

ijkm i'j'k'm' 

xafa'j'a^a^ l a m/ ak'a j/ ai'\0) 
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= E ^(l)^(l)^(2)^(2)^(3)^(3)^(4)V m (4) 

ijkm 

-6 E ^(l)^(l)^(2)^ fc (2)^ fc (3)^(3)^ ro (4)^ m (4) 

ifcj'm 

+8 2 ^(l)^(l)^(2)^ m (2)^(3)^(3)V m (4)^(4) (3) 

ikjm 

-6 2 ^(l)^ m (l)^(2)^(2)^(3)^(3)^ m (4)^(4) 

ikjm 

+3 2 ^(l)^(l)^(2)^(2)^(3)^(3)^ m (4)^(4). 

ikjm 

Let us mention that the space coordinates in Eqs.(2) and (3) coincide for all spinors 
in the ZRL and that the subscripts k in the last three lines of Eq. (2) and k,m in 
the last four lines of Eq. (3) run over single-particle occupied states. Let us emphasize also 
that the sum of coefficients of the exchange (Fock) terms and of the direct (Hartree) term 
in Eqs. (2) and (3) is equal to zero, as it should be. 

We define the scalar (p sq ), vector (p vq ) and tensor (pr q ) densities for neutrons (q = n) 
or protons (q = p) in the usual way: 

Psq (l)= E ^(l)Vi(l), 

i(q states) 

p vq (l)= E ^(1) 7 °^(1), (4) 

i(q states) 

pr ff (i)= E E ^(1)^(1)^ = 1 E ^(i)7°7-^(i), 

k=l i(q states) i(q states) 

where the subscript i runs over all the occupied states of a nucleon of type q. Here, n is the 
unit vector along the radial direction. 

The scalar, vector and tensor total densities are p s = p sn + p sp , p v = p vn + p vp and 
Pt = PTn + Pt p , respectively. 

Using the Fierz transformation [16], we can write for a fermion system: 



a, 6 a,b 

-{$alsl^a){4>blslM + ^(iv^)^^)]- ( 5 ) 

The quantities J2a,bA^^b)(Mc)(^c^ a ) and T,aMA 1 i J ^b)(Mc)(^c^d)(^d^a) can also 
be written in a similar way, although they involve much more terms and will not be given 
here. 
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In order to write down the quantities < 0|(?/>-?/ ; ) 3 |0 > an d < 0|(^>) 4 |0 > i n terms of 
the nucleon densities defined in Eqs. (4), we take into account the following relation: 



i k=l 

From this relation, one gets two more useful equations 



(6) 



and 



i i 

Y^A^^iYl^^i = Pt [(! - <W(1 - 8 K0 )n u n K - 5 m 5 k q\ 

i i 

Finally, from this last equation one obtains 



(7) 
(8) 



(9) 



Having in mind that the space coordinates in Eqs. (2) and (3) coincide for all spinors 
in the ZRL, we get from Eq. (2) for the cubic term 



(0|(^) 3 |0) 



= P\ 
1 

+ 8 



rsn * rsp * rvn * rvp 



(pTn + PTp) 



Pin + Pip + 3(PsnPln + Pspplp) ~ ^PsnPTn + PsppTp) 



and from Eq. (3) for the quartic term 



(10) 



(0|(^) 4 |0) 



Ps - nPs 



Psn Psp Pvn Pvp 



(PTn + PTp) 



-2 Ps 



32 



Pin + Pip + HPsnPln + Pspplp) ~ 3(psnpTn + PspP T 

4 , 4 , r>( 2 2,2 2 \ «/2 2 , 2 2 
Psn ' Psp ' V\PsnPvn ' PspPvp) G \PsnPTn ' PspPTp 



+Pvn + Pvp - 2 (PvnP T n + PvpPTp) + P T n + P T p 
I if 4- ff - if - ff - if, - ,•■-]- f 

16 



3 / . 

(p; 



2 

sn 



2 
Psp 



2,2 2 2 

Pvn ' Pvp PTn PTp 



Let us mention that the contributions given by Eqs. (10) and (11) include the densities 
p s , p v and pr but they do not include the contributions of pseudoscalar pp- and axial vector 
p^-densities for parity reasons. 
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From Eqs. (10) and (11), it can be easily seen that the exchange (Fock) contributions 
to the energy of the system corresponding to both the cubic and quartic terms are very 
essential. Actually, only the p^-term in Eq. (10) and the p^-term in Eq. (11) originate from 
the direct (Hartree) contributions, all the other terms in Eqs. (10) and (11) arise from the 
Fock contributions, and they have a strong isovector structure (ISVS) 1 . 

B.2 Contributions of the NLSI terms to the nucleon self-energy 

The contribution of the NLSI terms of the third (3) and forth (4) order to the self- 
energy of a nucleon of type q [S^ 3,4 -*] can be extracted from the following equation: 

E< 3 '% = ^U S T = [Eg' 4 > + Sg'V + i^fl°l ■ n]W (12) 

From Eq. (12), it is seen that the NLSI terms Ugi give a contribution to the total nucleon 
self-energy of the same structure as the linear components of the Lagrangian. Thus, each 
component of the total self-energy can be written as 

y-i s^linear i y(3,4) /i q\ 

^iq — ^iq ^iq ' 

where % — s, 0, T specifies the scalar, time component of the vector, and tensor self-energies. 
One can look in Refs. [3,14] for further details related to E'* near . 

Taking into account that 

5 5 5 

&f ps = ^' Jf pv = 7 °^' and &f pT = i7 ° 7 ' ^ 
one obtains for the contribution of the cubic terms of U S i to the self-energy components 



^sq = 3pl ~ ^ {PL + P% + P 2 vn + pl p ~ pTn ~ pip) ~ g (^PsPsq ~ p\ ~ pl q + pTq) , (15) 



S 3 g = \pvq{.Psq - 2p 8 ), (16) 



S T«J = ^PTq^Ps -Psq)- (17) 

For the contribution of the quartic terms of U si to the self-energy components, one has 



1 We consider that a quantity has ISVS if it depends on the difference p n — p p (or p sn — p sp ). 



= + 4 (Psq - A Ps) (pL + Pip + PL + P^p - PTn ~ PTp 

13 3 

+ ^{pln + Psp) + -^{Psnpln + PspP^) ~ g^W^n + PspPTp) (18) 

- 3 PsP s </ + ^P s (Ps 9 + P^ g - Pt 9 ) - gPs^Psq + 3 P^ - 3 Pt 9 )> 



T (4) - 3« 

^Og ~~ °Pvq 



^Tq — 3PT 9 



PsPsq -Pi- \(3plq + Plq ~ Pig) + \ (pL + Psp + PL + pip ~ PTn ~ Pt V ) 



PsPsq + ~(3p% + P^ 9 - Pig) ~ ^(PL + P^p + PL + P^ p - Pin - Pip) 



,(19) 



.(20) 



From the previous publications [3,14], it is clear that the isovector structure of the 
Hartree-Fock solutions is strongly determined by the contribution of the ir- and p-mesons. 
The present results, show that the isovector-independent NLSI also make essential contribu- 
tions to the ISVS of the energy and self-energies of the system. 

From Eqs. (12,13), it follows that the quantities Sg' 4) , S^' 4) and S^ 4) enter the RHF 
Dirac equation in the same manner as the self-energies Yi sq , S 0g and Er g , produced by linear 
interactions, do. However, the S^'^-components involve a strong density dependence, and, 
what is more surprising, a strong ISVS. The densities p sq , p vq and p Tq can be calculated in 
the following way 



M0= ^ E [G 2 {r)-FKr) 

i{q states) 

M0= ^ E [G 2 (0 + i?(r)] , (21) 

i(q> states) 

PT,(r)= ^ E ^(r)-F,(r), 

i(g states) 

where Gi{r)/r and Fj(r)/r are the radial functions of the upper and lower components of 
the nucleon Dirac spinor. Let us mention that Eqs. (10, 11, 15-17) are valid for NM and 
finite nuclei. 
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3 Numerical results 



A. Nuclear Matter 

We shall start the discussion of our results for symmetric NM, where px = 0, p SjP = p s , n and 
Pv,p — Pv,n- As is mentioned above, in the present paper we follow closely the Lagrangian 
of Ref. [14] corresponding to the set HFSI and, in our calculations, we have fixed the same 
parameters as in this reference. Thus, the bare nucleon mass and the 7r, oo and p meson 
masses have been taken to be equal to their empirical values: M = 939 MeV, m n = 138 
MeV, = 783 MeV, and m p = 770 MeV. As for the f n and g p coupling constants, we have 
chosen the experimental values fl/An = 0.08 and g 2 p /An = 0.55. The ratio f p /g p = 3.7 is 
taken in accordance with the vector dominance model. Then, we are left with the following 
five free parameters: m a , g a , g^, b and c, which are to be adjusted to reproduce some 
observables for NM and finite nuclei in a similar manner as in Ref. [14]. 

In Fig. 1, we present the result of the RHF calculation of the energy per particle (E/A) 
for NM as a function of p v for three cases: 1) The curve HFSI corresponds to the results 
obtained by the method suggested in Ref. [14] to take into account NLSI. 2) The curve 
ZRL*[b,c(HFSI)] corresponds to the results obtained with the Lagrangian and the method 
proposed in this paper, taking for the parameters m a , b and c the same values as in the HFSI 
approximation [14], g a and g u being chosen to get the saturation point at the same density, 
Po = 0.14 fm -3 , as in the HFSI set. We obtain in this case the compressibility modulus 
K = 319 MeV. 3) The curve ZRL[g a , g^, b,c(HFSI)] corresponds to the results obtained 
for the same set of parameters as in the HFSI set of Ref. [14], the NLSI being treated in 
the ZRL by the method suggested in the present paper (the saturation point is achieved at 
~ 0.113 fm -3 in this case) 2 . To check the reliability of the ZRL approximation, we have 
carried out calculations of the E/A vs p^-curves in the Hartree approach also. The results 
of the calculations are presented in Fig. 2. The ZRL approximation introduces appreciable 
differences with respect to the exact calculation for densities around the saturation value. 

Considering g a , g^, b and c as fitting parameters, it is easy to determine a set with a 

compressibility modulus K equal to a given reasonable phenomenological value. In Fig. 3, 

the energy per particle E/A as a function of p v is represented for the model considered in 

the present paper (ZRL) for three sets of parameters. In all sets the scalar mass and one 

coupling parameter are chosen with the same values as in the HFSI set of Ref. [14] , whereas 

the other three coupling parameters are chosen to get the same values of po an d E(p )/A in 

NM as in the HFSI set (i.e. p =0.14 fm~ 3 , E/A = -15.75 MeV) and 7^=250 or 275 MeV: 

2 In the limit of small densities, where the contribution of the self-interaction terms becomes negligible, 
the HFSI- and ZRL-curves coincide. 
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1) ZRL cl set: c is fixed to the value (c = —0.01461) obtained in the HFSI set, while 

— 2 2 

b = —0.011533, |^ = 3.0621, |^ = 7.849 generate the required saturation conditions of NM 
with K = 275 MeV. 

2) ZRL W set: 6 is fixed to the value (6 = -0.006718) obtained in the HFSI set, while 
c = 0.012475, f| = 3.3039, = 8.9346 generate the required saturation point of NM with 
K = 275 MeV. 

3) ZRL 62 set: b is fixed to the value (b = -0.006718) obtained in the HFSI set, while 

2 2 

c = 0.042561, |£ = 3.1715, |£ = 8.4777 generate the required saturation point of NM with 
K = 250 MeV. 

One could also try to get K = 250 MeV keeping the value of c as in the HFSI set, 
however, we could not find solutions in this case. From Fig. 3, one can see that appreciable 
differences between the three models appear only at densities larger than the saturation one. 
Thus, a similar description of finite nuclei can be expected for these three sets. 

Our next step is to find an adequate parameterization for NM and finite nuclei. To 
do that, we can follow the procedure used in Ref. [14]. In the present paper, the values of 
g a , b, c are determined by reproducing the saturation conditions of symmetric NM for po, 
E/A and a reasonable value of the compressibility modulus. In this work we have fitted 
two sets, ZRL1 and ZRL2, corresponding to the compressibility modulus 250 MeV and 275 
MeV, respectively. Calculations for finite nuclei put some extra constraints on the values of 
rricr and {m a is adjusted to get the experimental r.m.s. charge radius of the 16 and the 
value of g u is chosen to get reasonable values of spin-orbit splittings). Alternatively, one can 
use finite nuclei data to fit directly the free parameters of the model. 

The values of the parameters chosen in the present paper for the ZRL1 and ZRL2 sets 
and for the HFSI set of Ref. [14], as well as some calculated NM properties, are given in 
Table I. One can appreciate that the scalar meson mass m a needed to get a good description 
of the surface nuclear properties is considerably larger in the ZRL1 and ZRL2 sets than 
in the HFSI one. The parameters b, c are quite different in the sets ZRL1 and ZRL2 as 
compared with those of the set HFSI. This is due to the strong Fock terms contribution of 
the NLSI in sets ZRL1 and ZRL2. 

Fig. 4 shows for the ZRL1 set the values of E/A as a function of p v for different values 
of the asymmetry parameter 5 defined by the following equation: 



5 = 



Pv,n 



Pv,p 



(22) 



Pv,n 



~r Pv,p 
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It can be appreciated in this figure a strong shift of the equilibrium density po(5) 
towards smaller values as 5 increases and also that the pure neutron matter (5 = 1) appears 
to be unbound. The difference between the E/A values corresponding to the symmetry 
parameters 5 = and 5 = 0.2 remains almost constant for p v > 0.4 fm -3 . This fact is 
reflected in a strong change in the trend of the symmetry energy parameter 04 as a function 
of the density for p v > 0.4 fm -3 . This can be seen clearly in Fig. 5, where we have represented 
the CZ4 parameter in the ZRL1 set as a function of the density. The strong deviation of 
from the linear dependence at very high densities is related to the highly increasing role 
of the NLSI terms, especially the quartic ones, in this region. Anyway, this set generates 
reasonable solutions up to densities larger than the HFSI set does. In this latter case, to 
get a value of K — 250 MeV, one needs a contribution of the NLSI terms producing a rapid 
decreasing of m* with the density for p v > p [14]. Thus, m* becomes almost zero for 
densities larger than p v ~ 0.22 fm -3 . That is why the quantities E/A and a 4 have only been 
represented up to this density in Figs. 1 and 5. 

For comparison, we have also represented in Fig. 5 with the NLHF label the results for 

04 obtained in Ref. [10], within an approximate procedure to include exchange contributions 
of the NLSI terms of scalar type. One can see from this figure that the results corresponding 
to the NLHF approximation lie significantly below our results. However, the NLHF approx- 
imation does not include isovector mesons. Of course, the inclusion of these mesons would 
significantly increase the value of 04. 

We have mentioned above that the NLSI terms are explicitly isovector-independent 
[see Eq. (1)]. However, their contributions to the nucleon self-energies manifest a strong 
isovector structure due to the exchange part. 

To illustrate this point, we show in Figs. 6-9 the density dependence of the quantities 
Ej,£, S^, Sq^, Sq£, respectively. The first two quantities correspond to the exchange 
(Fock) contributions of the cubic and quartic terms, respectively, of the NLSI part of the 
Lagrangian to the scalar nucleon self-energy, whereas the last two quantities represent the 
exchange contributions to the (time component of the) vector nucleon self-energy. We have 
considered three values of the asymmetry parameter: 5 =0, 0.5, 1. Notice that for the case 

5 = 1, the proton self-energies correspond to a proton moving in pure neutron matter and 
that, in this case, S^ 4 '* are equal to zero. 

From Figs. 6-9, it is seen that the exchange (Fock) self-energies S-^ 4 -* (i = s, 0) 
represent a very important contribution to the respective total self-energies even at normal 
densities of NM (especially that of £^ 4 ' ) . This contribution is appreciably smaller for the 
ZRL2 set with a larger value of the compressibility modulus, due to the smaller values of 
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the parameters b and c. It is also seen that (i — s, 0) are essentially dependent on the 

neutron excess in the system (due to the Hartree-Fock framework used). 

B. Finite nuclei 

In order to make a more complete analysis of the properties of our nonlinear RHF 
model, we have carried out calculations for finite nuclei with the ZRL1 and the ZRL2 sets of 
Table I. Actually, as we have explained above, we have taken into account, besides the satu- 
ration NM properties, the experimental values of the binding energies, spin-orbit splittings 
and r.m.s. charge radii of finite nuclei in choosing the parameters of the ZRL1 and ZRL2 
sets. We remind that the ZRL1 and ZRL2 sets contain five fitting parameters, as explained 
in subsection A dedicated to NM. 

In Tables II and III, we present the RHF results of our ZRL1 set 3 , containing the 
exact exchange contribution of the NLSI, given by Usi in Eq. (1), for the ground state 
properties for the five doubly-magic nuclei: 16 0, 40 Ca, 48 Ca, 90 Zr, 208 Pb. For comparison, 
we also present the results of our HFSI set of Ref. [14] (where the exchange terms of NLSI 
are treated in an approximate way) and experimental values. The CM. values in Table II 
indicate an estimation for the center-of-mass correction to the total energy, which is not 
included in the ZRL1 and HFSI results. 

From Tables II and III, one can see that the ZRL1 set allows a quite good description 
of binding energies, spin-orbit splittings and r.m.s. charge radii for spherical nuclei. The 
results are comparable to those of the HFSI set, although, as we have already indicated 
above, the ZRL1 set has reasonable solutions for NM up to much higher densities than the 
HFSI one has. As is common in the RHF approaches containing the exchange of pions [9], 
the ZRL1 set also predicts a strong reduction of the spin-orbit splitting as going from the 
40 Ca nucleus to the 48 Ca one. We explained in Ref. [9] that this fact is a consequence of 
the small value of the pion mass. 

Figs. 10-14 show the calculated charge distributions for the indicated nuclei and the 
corresponding experimental ones for comparison. It is remarkable that there is a very good 
agreement between the theoretical results and the experimental data, specially for the 16 
and 40 Ca nuclei. 

3 Thc results for finite nuclei with the ZRL2 set are almost identical to those obtained with the ZRL1 set. 
This is due to the fact that the properties of NM obtained with both paramctrizations are very similar in 
the range of densities relevant in finite nuclei. Only for larger densities, the properties of these two sets are 
significantly different. Notice that the smaller values of the b and c coefficients in the ZRL2 set than in the 
ZRLf one is related to the larger value of the compressibility modulus in the former case than in the second 
one. 



12 



4 Conclusions 



In this work, we have considered a Lagrangian including, beside the exchange of a, u, it and 
p mesons between nucleons, NLSI of zero range, associated with three- and four-body forces, 
in the framework of the RHF approximation. The Fierz transformation allows us to express 
the exchange (Fock) terms of the NLSI explicitly via the direct (Hartree) terms in an exact 
way. The model has been applied to NM and finite nuclei. 

The NLSI (4>vp) 3 and (t/>-?/>) 4 considered here are constructed from the scalar densi- 
ties. However, the relativistic Hartree-Fock procedure generates, in this case, the nucleon 
self-energies (in the respective single-particle Dirac equation) with the space-time trans- 
formation properties of a relativistic scalar, vector and tensor (the pseudoscalar and axial 
components do not survive). It should be mentioned also that the NLSI generate strong 
density dependence of the respective nucleon self-energies £^' 4 \ £^' 4 \ and X^ 4 "*- m the 
RHF framework, the method suggested manifests also a strongly developed ISVS of the self- 
energies generated by isovector-independent NLSI, this point being one of the basic features 
of the RHF framework (notice that in the relativistic Hartree approximation there is no 
ISVS due to the NLSI, neither, in the usual form a 3 - and <r 4 nor in the ZRL). 

Let us mention that the same method, as suggested in the present paper, can be utilized 
to treat the self-interactions of the o/-type (u being the vector-isoscalar meson field) and 
also associated with four-body forces [15,17]. In this case, in the ZRL, the nonlinearities 
are constructed by mean of the relativistic vector densities ip^ip and contain eight fermion 
operators. Self-interactions of the cj 4 -type will introduce an additional fitting parameter. We 
shall consider this case in more detail in a forthcoming paper. Actually, NLSI of different 
types (<ru; 2 , a 2 u 2 , etc. [15]) can be treated in the ZRL in the same fashion as we did for the 
Us i term. 

The results given in sect. Ill show that our model with five parameters is flexible enough 
to generate NM properties, in the range of densities relevant for finite nuclei, compatible with 
data inferred from nuclear experiments. Although the present model allows to get reasonable 
results for NM until values of the density much higher than those of the HFSI set of Ref. [14], 
we cannot expect that our results, obtained with a parametrization that describes adequately 
finite nuclei properties, could be considered realistic for high densities, for example, larger 
than 2p - 

The calculations for doubly- magic nuclei have been also carried out within the ZRL1 
set. Our results show a rather good agreement between theory and experiment for binding 
energies, spin-orbit splittings, r.m.s. charge radii and charge densities, especially, in light 
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and mid-weight nuclei. 
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Table 1. Adjusted parameters and some properties of symmetric NM for the parameter sets ZRL1, 
ZRL2, corresponding to the zero range limit (ZRL) considered here, and HFSI from Ref. [14]. The 
value of the effective nucleon mass M* is defined in the conventional way and is given at the Fermi 
surface. 04 is the symmetry energy parameter. The equilibrium density po is given in fm~ 3 , the 
mass of the scalar meson m a , the energy per particle at the equilibrium E/A, the compressibility 
modulus K and the symmetry energy parameter 04 are given in MeV. 

Set gl/4-K gl/4n m a 6 * 10 3 c * 10 3 p E/A K M*/M a 4 

ZRL1 5.5743 11.732 497.8 2.9646 51.00 0.155 -16.39 250 0.58 35.0 

ZRL2 5.5121 11.565 496.8 2.4591 46.50 0.154 -16.37 275 0.58 35.0 

HFSI 4.005 10.4 412.0 -6.718 -14.61 0.140 -15.75 250 0.61 35.0 



Table 2. Comparison of the results of the present calculation for the ZRL1 set in finite nuclei 
with the corresponding results of the HFSI set of Ref. [14] and the experimental ones. The total 
binding energy per particle E/A, the non-relativistic center-of-mass correction to E/A, and the 
proton spin-orbit splitting Als for the shells lp of the 16 O nucleus and 2d of the 40 Ca and 48 Ca 
nuclei are given in MeV (the experimental values of Als in the Ca isotopes are not very well 
established [9]). The r.m.s. charge radii r c are given in fm. 

16 A0 Ca 48 Ca 90 Zr 208p b 

Set -E/A r c A LS -E/A r c A LS -E/A r c A LS -E/A r c -E/A r c 

ZRL1 7.37 2.71 6.3 8.33 3.44 7.1 8.51 3.49 2.6 8.67 4.25 7.85 5.49 

HFSI 7.43 2.73 6.4 8.33 3.48 7.05 8.45 3.48 3.27 8.58 4.26 7.78 5.52 

EXP 7.98 2.73 6.3 8.55 3.48 6-7.6 8.67 3.47 5 8.71 4.27 7.87 5.50 

CM. 0.61 0.20 0.18 0.08 0.02 
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Table 3. The single particle energies (in MeV) for protons and neutrons in the 16 0, 40 Ca, and 
48 Ca nuclei. For each state, the first and the second rows correspond to sets ZRL1 of the present 
paper and HFSI of Ref. [14], respectively. References to the experimental data are given in Refs. 
[9] and [14]. 

ie O 40 Ca 48 Ca 



protons 



neutrons 



protons 



neutrons 



protons 



neutrons 



State calc. exp. calc. exp. cole. exp. calc. exp. calc. exp. calc. exp. 



lsi/2 39.51 40±8 43.89 47 



38.96 



43.15 



48.71 50±11 57.22 
48.64 56.83 



54.22 55±9 59.01 



55.3 



58.8 



lp 3/2 19.28 18.4 23.36 21.8 33.39 34±6 41.53 
18.80 22.73 32.2 40.08 



39.33 35±7 43.08 
39.7 41.95 



lp 1/2 12.98 12.1 16.95 15.7 29.43 34±6 37.47 
12.40 16.22 27.8 35.56 



38.24 35±7 41.24 
38.0 39.60 



Id 



5/2 



17.22 14.3-16 24.98 21.9 23.04 20.5 26.20 16 
16.2 23.77 23.2 25.0 



2s 



1/2 



8.64 10.9 16.36 18.2 15.30 15.8 18.42 12.4 
9.83 17.23 15.8 18.5 



Id 



3/2 



10.10 8.3 17.71 15.6 20.45 15.5 21.74 12.4 
9.19 16.55 19.9 20.2 



If 



7/2 



9.98 9.9 
9.21 
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Fig. 1: The E/A values for symmetric NM in the RHF approach in three cases: The HFSI (dashed) 
curve corresponds to the results obtained with the HFSI set of Ref. [14]; the ZRLf^o-, <7 w ,b,c(HFSI)] 
(solid) curve corresponds to the results obtained by the method suggested in the present paper, 
with the same set of parameters as in Ref. [14]; and the ZRL*[b,c(HFSI)] (dash-dotted) curve 
corresponds to the results with the parameters g a and g u chosen to get saturation at the same 
point as in the HFSI set. 
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Fig. 2. The E/A values for symmetric NM in the Hartree approach in two cases: The NL3 
(dashed) curve represents the results obtained with the NL3 parametrization of Ref. [20] and the 
ZRL[H:NL3] (solid) curve are the results obtained with the NL3 parametrization in the ZRL. 
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Fig. 3. The E/A values for symmetric NM, with the Lagrangian considered in sect. I (which 
contains the zero range NLSI given in Eq. (1)) for different choices of the parameters fitted to the 
same values of NM po and E{po)/A as in the set HFSI of Ref. [14] (see the text), and generating 
two values of the compressibility modulus: The ZRL cl (solid) curve, with c < taken from the 
HFSI set, corresponds to K = 275 MeV; the ZRL 61 (dashed) curve, with b taken from the HFSI 
set, corresponds to K = 275 MeV; and the ZRL 62 (dash-dotted) curve, with b as for the ZRL 61 
case, corresponds to K = 250 MeV. 
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Fig. 4. The E/A values in asymmetric NM, for the ZRL1 set, with different values of the asymmetry 
parameter 5. 
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marcoss 1-Sep-2003 10:55 

Fig. 5. The symmetry energy parameter 04 as a function of the NM density for sets ZRL1, HFSI 
[14] and NLHF [10]. 
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Fig. 6. The density dependence of the exchange (Fock) contribution of the (tpip) 3 term to the 
scalar nucleon self-energy for different values of the asymmetry parameter 5. 
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Fig. 7. The same as in Fig. 6 but for the density dependence of the exchange (Fock) contribution 
of the (^V) 4 term. 
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Fig. 8. The density dependence of the exchange (Fock) contribution of the (tptp) 3 term to the 
(time component of) the vector self-energy. 
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Fig. 9. The same as in Fig. 8 but for the density dependence of the exchange (Fock) contribution 
4b of the (V^) 4 term. 
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Fig. 10. The charge distribution (dashed line) for the 16 nucleus, with the ZRL1 set. The 
experimental charge distribution (solid line [21]) is also shown for comparison. 
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